On the Incommensurate Phase in Modulated Heisenberg Chains 



OO 



F. Schonfeld, G. Bouzerar, G.S. Uhrig and E. Miiller-Hartmann 

Institut fur Theoretische Physik, Universitdt zu Koln, 
Ziilpicher Str. 77, Koln 50937, Germany. 
(February 1, 2008) 

Using the density matrix renormalization group method (DMRG) we calculate the magnetization 
of frustrated S = 5 Heisenberg chains for various modulation patterns of the nearest neighbour 
coupling: commensurate, incommensurate with sinusoidal modulation and incommensurate with 
solitonic modulation. We focus on the order of the phase transition from the commensurate dimer- 
ized phase (D) to the incommensurate phase (I). It is shown that the order of the phase transition 
depends sensitively on the model. For the solitonic model in particular, a fc-dependent elastic energy 
modifies the order of the transition. Furthermore, we calculate gaps in the incommensurate phase 
in adiabatic approximation. 
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I. INTRODUCTION 

Low dimensional magnetic systems have attracted con- 
siderable attention in recent years. Various theoreti- 
cal and experimental efforts have been made to under- 
stand the fascinating low energy physics of quasi-one di- 
mensional gapped spin systems, such as spin-Peierls sys- 
tems ( CuGe0 3 or NaV 2 A @||), Haldane systems 
(e.g. Ni(C 3 HioN2)2N)2C104 §) and spin ladders (e.g. 
SrCu 2 3 § orCu 2 (C 5 Hi 2 N 2 )2Cl4 §). Even if questions 
still remain open, many of the experimentally observed 
features can be already understood within the framework 
of one dimensional Heisenberg chains with various cou- 
plings (this includes also spin ladders [Q). 

Some of these systems exhibit interesting features in 
external magnetic fields, for instance, a transition from 
a commensurate to an incommensurate phase. At this 
transition weak hysteresis effects are observed in CuGc0 3 
at low temperatures |^|| and Kiryukhin et al. found a 
small jump in the incommensurability measured by X-ray 
scattering . These features are characteristic for a 

first order phase transition. From the theoretical point of 
view, no consensus has been reached so far on the order of 
the transition. The phase transition was predicted to be 
of first order by Cross |L3]]. Bhattacharjee et al. obtained 
the same conclusion using a phenomenological Landau 
expansion [|l4j. But mean-field calculations of Fujita and 
Machida for a renormalized XY-model display a second 
order phase transition Jl[| while Buzdin et al. [[l6) find a 
second order phase transition only at T — and a first 
order one for T > using essentially the same model 
as Fujita/Machida. Horovitz underlines the importance 
of the correct treatment of cutoffs when passing to the 
continuum limit jl7],[l8| ■ 

In this paper, we propose to clarify which parameters 
influence the properties of the I phase with the help of the 
DMRG method for finite systems. In section II we cal- 
culate magnetizations for different types of modulations 
and show that the order of the D-I phase transition is 
model dependent. In the I phase we calculate the magne- 
tization dependence of the two gaps A + /_ corresponding 



to the increase/ decrease of the z component of the total 
spin by unity Jl9|. 

For all calculations we have chosen parameter sets 
which are convenient for the numerical calculations, i.e. 
displaying small finite size effects. Computational as- 
pects are given in section III. In section IV we summarize 
the results. 



II. MAGNETIZATION 

In the adiabatic approximation for the phonons the 
modulation of the exchange couplings can be described 
by parameters Si which are linked to the lattice distor- 
tion. Thus the Hamiltonian includes an elastic energy 
which is a positive quadratic form of the {Si}. In a first 
step we take the elastic energy to be dispersionless, i.e. 
diagonal in real space 
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where a denotes the relative frustration and S z is the z 
component of the total spin of the L-site chain. The last 
two terms in (Q) are the Zeeman energy and the elastic 
energy associated to the lattice distortion. 



A. Fixed Modulations 

As a starting point let us consider the simple case 
where the lattice distortion is kept frozen as in the D 
phase even in the presence of a magnetic field. 



Si = (-iys 
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The amplitude S is treated as a fixed parameter. The 
constant elastic energy is not taken into account for the 
moment. Chitra and Giamarchi |^0| calculated the mag- 
netization of frustrated or dimerized spin chains in a 
magnetic field using bosonization techniques. Within 
this continuum- limit approach the frustration and dimer- 
ization cannot be treated simultaneously (double sine- 
Gordon model). For a < a c and 6 > the contribu- 
tion of the frustration is assumed to be irrelevant and 
the model reduces to an integrable sine-Gordon model. 
However, recently it has been shown by Bouzerar et al. 
that this is not the case |2l| . The magnetization increases 
just above the lower critical field H c like m oc \/H — H c 
p0| , p2| . The same power law is found for a > a c and 
5 = 0. In particular, the transition to finite magnetiza- 
tion is of second order in both cases. With the DMRG 
method, we find this square root behavior in presence of 
both dimcrization and frustration. But as shown recently 
by Tonegawa et al. by means of exact diagonalization 
p3[ an additional remarkable difference appears in some 
parameter range of dimerization and frustration, i.e. a 
plateau at m = 1/4, see Fig. |l|. 



a constant contribution at given amplitude which will be 
dropped for the following consideration. 

Using the Jordan- Wigner transformation the applied 
magnetic field corresponds to a shift of the chemical po- 
tential. For the XY— model with a finite magnetization 
m = S z /L, it is straightforward to show that an infinites- 
imal spin-lattice coupling leads to an instability at mo- 
mentum q = 2kp = 7r(l + 2m). In the case of the Heisen- 
berg model, this relation is expected to hold true as well 
Jl3|,^9 20 2j],|2(|. We have confirmed numerically for var- 



ious sets of parameters <5, a and various system sizes that 
the energy is minimum at q = tt(1 + 2m) for a given m. 
The ground state energy per site for H = as a func- 
tion of q is plotted in Fig.||, for various magnetizations 
m at fixed S = 0.1 and a = 0.35. The positions of the 
cusps correspond exactly to q = tt(1 + 2m). The fact 
that one observes cusps and not smooth quadratic min- 
ima is linked to the divergence of the susceptibility at 
q = 7r(l + 2m), i.e. an instability even for infinitesimal 
coupling. This is the generic behavior, independent of the 
parameters a and 6, and confirms the relation between 
the wave vector and the magnetization. 




o 
LU 



-0.30 



-0.32 



-0.34 



-0.36 



-0.38 



-0.40 




FIG. 1. Magnetization as a function of the magnetic field 
for 8 = 0.3, q = 0.1 (filled circles) and S = 0.5, a = 0.2 (open 
squares) for a 60 site chain. 



FIG. 2. Ground state energy per site of an 80 site chain as 
a function of the wave vector q for various magnetizations m, 
a = 0.35 and 8 = 0.1. 



In the case of finite magnetization, it is known that 
the commensurate dimerization pattern (fl) is not appro- 
priate for describing spin-Peierls systems. For instance 
X-ray measurements on CuGeOs clearly show that the 
structure of the lattice distortion becomes incommen- 
surate under a sufficiently large magnetic field jio] |l2| . 
Thus a more appropriate choice for the modulation is, 

= (Jcos^) (ii) 

as it was suggested in ]l9| , p4[ . To begin with, q is con- 
sidered as a free parameter which is fixed by minimizing 
the total (free) energy. Note that for 5 / 7r the elastic 
energy is q independent for the ansatz (y) yielding only 



Henceforth, we fix q = 7r(l + 2m) and investigate 
the magnetization as a function of the applied field. 
We find that the incommensurate exchange coupling ([^) 
has a rather strong effect on the magnetization lead- 
ing to a first order phase transition. To elucidate this 
we present the magnetic ground state energy per site 
E(m) = ((ffchain) + E c \ &st )/L as a function of the mag- 
netization in FigJ^. Results for several chain lengths are 
included to show the absence of finite size effects. 

The salient feature of E(m) for sinusoidal modulation 
is the discontinuous jump at m = 0. To understand this 
jump it is helpful to look at the averaged squared distor- 
tion J^i $i which takes the value 5 2 at q = ir and S 2 /2 



2 



otherwise. We see that in the D phase all Si are max- 
imally distorted whereas in the sinusoidally modulated 
phase there are also large regions with weaker distortion. 
So neither the elastic energy is not continuous in the limit 
q — ► 7r nor is the magnetic energy since it reacts to the 
distortions. 
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FIG. 3. Open symbols (left scale): ground state energies 
E(m) — E(0) as a function of the magnetization for the sinu- 
soidal modulation (0) (8 = 0.2, a = 0.35) for chains of 100 
(circles) 80 (squares) and 60 (diamonds) sites. To highlight 
the discontinuity at m = a cubic fit for m > is depicted 
with a solid line. The inset shows an enlargement and the 
tangent for m — 0.05 as described in the text. 
Filled circles (right scale): ground state energies E(m) — E(0) 
for the adaptive modulation from (jli|), Ko = 1.7 (S ~ 0.2 in 
the D phase) and a = 0.35. The dot-dashed line is just a 
guide to the eye. 

To deduce the dependence m{H) from FigJ| we have 
to resort to Maxwell's construction, i.e. we compute the 
convex hull. The magnetic field defines the slope g^H — 
dE/dm of the tangent which touches the convex hull at 
the value m (Legendre transformation). So one obtains 
m(H). The jump in E(m) leads to a first order transition 
with a jump in m(H). The resulting m(H) deduced from 
Fig.|| is depicted in Fig.f|. 

Calculating the corresponding local magnetizations 
|jl9f one finds that there is a large alternating local mag- 
netization close to each zero of the modulation. Summing 
the local magnetizations around each zero one finds a 
contribution of S z = 1/2, i.e. of one spinon. 
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FIG. 4. Magnetization as a function of the applied mag- 
netic field for a = 0.35 and 5 = 0.2 of a 80 site chain as 
deduced from the open squares in Fig.^. The inset shows 
an enlargement near H c for 100 (circles) 80 (squares) and 60 
(diamonds) site chains. 



B. Adaptive Modulations 

In the previous section we chose a sinusoidal modula- 
tion and found that q — 7f(l + 2m) minimizes the total 
energy. We now proceed in a more general way by mini- 
mizing the total energy including the elastic energy term 
with respect to all the parameters Si. In other words, 
we allow the lattice distortion to adapt to the spin sys- 
tem. Within our DMRG approach we follow the iter- 
ative procedure proposed by Feiguin et al. |2^] who ap- 
plied exact diagonalization and Monte-Carlo-Simulations 
to a slightly different model. The dimerization ampli- 
tudes Si are calculated self-consistently by minimizing 
(-Hchaitt) + ^oiast which leads to 



J{8i-Si+i)+K 6i 



J 



5> 
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where the last term ensures that the 5, satisfy the con- 
straint Yli=i $i = 0- Following j28| this equation is used 
to improve iteratively the local distortions Si. The expec- 
tation values are taken with respect to the ground state 
of the previous iteration. Starting from the sinusoidal 
modulation (0) we find that four to five iterations are 
enough to achieve a stable pattern that does not change 
significantly on further iterations as shown in Figj^. 



3 



0.2 



-0.2 



» sinusoidal 



4. iteration 




10 



20 



30 



40 



FIG. 5. Incommensurate modulation: local distortions &i 
versus site index i for the first four iterations starting from 
a sinusoidal modulation (filled circles); 100 sites, a = 0.35, 
Ko — 3.3, and S z = 1. 

The envelope of the final modulation can be fitted by a 
product of complete Jacobi elliptic functions of modulus 
k as predicted analytically fll^ , |l(||l5| ] . For very low mag- 
netization, i.e. low concentration of solitons, the vicinity 
of each zero resembles a tanh [^7j . 

Within the self-consistent approach we calculate E(m) 
per site as plotted in Fig.|| (filled circles). E(m) is convex 
but in contrast to the curves with fixed sinusoidal modu- 
lation it is continuous. The convexity will be shown more 
clearly below in a modified representation. Thus we have 
a continuous, second order transition from the D phase to 
the adaptive I phase. The corresponding magnetization 
m(H) is shown in Fig|j| (filled circles). 

The enormous steepness of the continuous magneti- 
zation is explained by the following argument. For non- 
interacting spinons which are far enough from each other, 
the energy per site E(m) — E(0) is proportional to the 
number of spinons and hence to the magnetization m 
(see also filled circles in Fig.|| for small m). The pro- 
portionality constant eo is the energy of a single spinon 
and determines the critical field 2eo = gn B H c , since two 
spinons are created by breaking one singlet. Because 
the spinons are exponentially localized (cf. Fig.||) two 
spinons at mutual distance I have additionally an expo- 
nential interaction w(l) = wq exp(— cl). Here c is a con- 
stant of the order of the inverse correlation length and wq 
is a proportionality constant which is positive for repul- 
sion and negative for attraction. The typical distance of 
the spinons is I = l/(2m) since each spinon carries spin 
S = 1/2. Hence for not too large values of m the total 
energy in an external magnetic field H equals 



E(m) - E(0) = g^J, B {H c - H)m + w a 2me ^ 



(2) 



By minimizing this expression for repulsion (wo > 0) one 
derives H(m) which increases exponentially slowly just 
above H c . This in turn leads to the drastic increase of m 
as depicted in Fig|| 
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FIG. 6. Filled circles: Magnetization for the adaptive mod- 
ulation (Ko = 1.7, a = 0.35) as deduced from the correspond- 
ing curve in FigM 

Open squares: Magnetization with a dispersive elastic energy 
as discussed at the end of the paragraph. K = 6Kq, Ko and 
a unchanged. 

To present the effects of soliton interaction more clearly 
we pass to an affine representation of the ground state 
energy E(m) by investigating 



E^(m) := E(m) - E(0) - gii B H c m 



(3) 



which would be constant zero if no interaction between 
the solitons existed. Note that E c g(m) is convex if and 
only if E(m) is. In Fig. [7] the generic resulting curves are 
shown with filled symbols (solid line) for the XY— model 
and the spin isotropic XYZ— model. The results for the 
XY— model are obtained for an infinite system without 
frustration by a continued fraction technique based on 
Green's functions [fisff . Using this method the simpler 
solvability of the XY— model allows to iterate up to 80 
times for an infinite chain with periodicities up to 120. 
These data are included as an additional check that no 
spurious effects due to finite size or insufficient iteration 
are investigated. 

The results in FigJ^ for a dispersionless elastic energy 
comply perfectly with exponentially repulsive solitons 
(||). There is no sign of a long range interaction oc l/l 
as postulated by Horovitz for finite cutoffs as they occur 
naturally in discrete lattice models |l7],[ll|. In particular, 
no attraction for dispersionless elastic energies are found 

A dispersionless elastic energy is of course a drastic 
simplification of the real phononic system. Cross already 
argued jl3| that a pinning in fc-space should influence 
the order of the transitions. We expect that the D— >I 
phase transition becomes first order if the elastic energy 
itself favors the distortion at k = n. This means that 
K(n) is minimum if the elastic energy can be expressed 
as i?ciast = | J2k K(k)\5k\ 2 - The argument compares for 
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FIG. 7. Affine representation of the ground state energy. 
The long-dashed lines indicate the convex hulls to the lower 
curves. XY — model: The upper solid curve shows E e g for a 
dispersionless elastic energy with Ko = 0.625. The lower solid 
curve shows -E e ff for a dispersive elastic energy (K = 6Ko) as 
discussed in the following section. Both curves are obtained 
via the continued fraction technique. The filled and open 
diamonds depict DMRG results for an 80 site chain. 
XYZ— model: DMRG results in the dispersionless case (filled 
circles) and for K = 6A"o (open squares) for Ko = 1.7 and 
a = 0.35. The dashed lines are guides to the eye. 

a given wave vector q close to tt the elastic energy K(q) 
of a sinusoidal modulation (£j|) with the one of an array of 
domain walls with the same periodicity 2n/q. Since the 
latter has also contributions of higher harmonics ±3<7, 
±5q, ±7q, ... its elastic energy is J2 n l a (2rn-i) 9 | 2 ^(<7) 
where the coefficients |a( 2n +i)g| 2 are symmetric about n. 
Thus the elastic energy is higher than the one for the 
sinusoidal modulation. By this mechanism higher har- 
monics are suppressed due to the elastic energy leading 
to a smoother and more sinusoidal modulation. If the 
convex curve for the adaptive modulation in Fig.^J is in- 
fluenced in a way to approach the discontinuous curve for 
sinusoidal modulation one must expect a region of con- 
cavity for small m. Hence the convex hull differs from the 
curve itself and a jump in the magnetization occurs: the 
transition is first order. Put differently, we expect that 
a dip in the elastic energy at the zone boundary leads to 
an attraction of the solitons. 

To investigate our hypothesis numerically we use 
K(k) — K + 2Kcos(k) with K = K - 2K kept fixed 
to refer to the same amplitudes in the D phase. This 
elastic energy corresponds in real space to 
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(4) 



where 8 is a vector with components Si and K is a cyclic 
tridiagonal Lx L symmetric matrix of coupling constants 



with diagonal elements K and off-diagonal elements K. 
Generic results for the energies E(m) in affine representa- 
tion are depicted with open symbols (solid line) in FigJ^. 
We find indeed a concavity for small values of m. This 
implies soliton attraction and a first order transition. 

Furthermore, we show in Fig.|| the resulting magneti- 
zation curves with and without dispersion of the elastic 
energy. The difference between the second order tran- 
sition for the elastic energy without dispersion and the 
first order transition with dispersion is clearly visible. 
Additionally, the critical field H c at which the transition 
occurs rises on inclusion of the dispersion. This complies 
also with the above consideration since the energy of a 
single soliton rises due to K + 2K cos(fc) > Kg except at 
k = tt for K > 0. 

Finally, in Fig.|| the modulation patterns with and 
without dispersion are compared. Indeed, the inclusion 
of K > makes the modulation softer and more sinu- 
soidal. In conclusion our numerical results convincingly 
corroborate our expectations for the effect of a dispersive 
elastic energy. 
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FIG. 8. Modulations for S z = 1 for the same parameters 
as in FigJ^ for the XYZ— model. 

Numerically, we are not able to decide whether an ar- 
bitrarily small K already yields a soliton attraction. For 
smaller values of K the minima in the affine represen- 
tation occur for smaller and smaller magnetization and 
they are more and more shallow. We expect that the 
soliton attraction exists down to arbitrarily small values 
of K but it may become irrelevant in practice due to the 
exponential smallness of the corresponding energies. 

We also investigated negative values of K. No qualita- 
tive change of the soliton interaction was found in com- 
parison to the dispersionless case. The iterative proce- 
dure, however, becomes quite unstable already for small 
negative values of K. 



C. Adiabatic Gaps 

So far we aimed at the average magnetization as a func- 
tion of the applied magnetic field. Another interesting 
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quantity which is accessible once E(m) can be computed 
are the adiabatic gaps. It is a so far unsettled question 
whether spin-Peierls systems have or have not gaps in the 
incommensurate phase. 

On the one hand, it seems clear that the incommensu- 
rate modulation pattern can be shifted along the chains 
without energy cost. This is certainly true in the conti- 
nuum description and thus most probable also for not too 
small correlation lengths. This quasi-continuous symme- 
try gives rise to quasi- Goldstone bosons called phasons 
which are gapless . They do not change the spin sec- 
tor and thus have AS Z = 0. The physics of phasons is 
beyond an adiabatic treatment of the lattice distortion 
since within an adiabatic treatment the distortion is as- 
sumed to be fixed. 

A different issue is the question whether the gaps A± 
corresponding to AS Z — ±1 are finite or not. Note that 
these gaps do not need to be equal since the spin rotation 
symmetry is broken for finite magnetization. From a non- 
adiabatic viewpoint one can still infer from the smooth- 
ness of the E(m) curves that there are no such gaps in 
the I phase since the modulation adapts always to the 
average magnetization. Applying, however, an operator 
like S + (k) or S~(k) p5| , |f9[ and asking for the accessi- 
ble excitation spectrum may lead to a different answer. 
These operators act only on the spin part of the ground 
state and leave the modulation unchanged. Thus it is not 
unreasonable to expect that the gapless excitations are 
not accessible if their access required a re-arrangement of 
the whole, in reality three-dimensional, modulation. The 
underlying question is whether the states S' ± (fc)|5 2 ) are 
orthogonal to \S Z ± 1) or not, if we denote by \S Z ) the 
ground state for the magnetization S z . 

Here we will investigate the simpler question whether 
in the strictly adiabatic framework the gaps A + are finite 
or not. Uhrig et al. were only able to compute A + + 
A_ since this quantity did not require the knowledge of 
the corresponding magnetic field. 

We define by E(m, H) :— E(m) — mg^H the ground 
state energy with self-consistently optimized modulation 
{Si}- By 



E ± (m,H) :-- 
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we denote the ground state energy with one additional 
spin flipped upward (+) or downward (-), respectively, 
but with the modulation {5i} belonging to S z = mL, 
not to S z = mL ± 1. This means that for E±(m, H) the 
modulation is not optimized for the given magnetization. 
This corresponds to the situation accessible by applica- 
tion of S + (k) or S~(k) without reaction of the lattice 
part. Then the gaps are defined by 



The gaps A + and A_ for a 100 site ring are displayed 
in Fig.|| for a — 0.35 and K = 2.38 corresponding to 
5 w 0.14 in the D phase. Finite size effects are not 
yet completely negligible, but the qualitative behavior 
is the one shown and is in agreement with previous self- 
consistent renormalized Hartree-Fock results fll9l. 




0.12 

FIG. 9. The energy gaps A + (filled squares) and A_(open 
squares) as a function of the magnetization for L = 100, 
a = 0.35 and K = 2.38 

Most importantly, we can show by Fig.^ that both gaps 
are indeed finite and of equal order of magnitude. It is 
interesting that apparently A+ is smaller at small mag- 
netization and A_ is smaller at larger magnetization. At 
least in the adiabatic approach, we can show that even 
the I phase has gaps. It would be interesting if any exper- 
imental evidence in favor of the existence of these gaps 
was found. 



III. COMPUTATIONAL DETAILS 



A±(m) = E±(m, H) - E(m, H) . 



(0) 



In our DMRG calculation [ |30[|3l| ] we apply periodic 
boundary conditions to minimize finite size effects. We 
keep 128 (64) states in the truncation procedure. To 
account for the incommensurate structure we use the fi- 
nite size algorithm |^,^l[. In the first steps where the 
system is iteratively increased we use the reflection of 
the left hand block to build up the superblock although 
the reflection symmetry is not given at this stage. This 
initial error is reduced either by supplementary sweeps 
through the system of the desired length or by interme- 
diate sweeps through the system the length of which is 
commensurable with the lattice modulation. We tested 
the accuracy of the DMRG results by comparing the low- 
est energies of the XY- Model with sinusoidal modulation 
in different S^-subsectors with exact results for a 60-site 
ring. Keeping 128 (64) states we find the typical error to 
be smaller than 10~ 6 (10~ 5 ) for S z = and 10~ 5 (10~ 4 ) 
in higher S'z-subsectors. For the selfconsistent calcula- 
tions in the case of adaptive modulations we used the 
sinusoidal modulation (|ii|) as a starting configuration for 
the curves presented. For larger dimerizations, however, 
it is more convenient to start with a step-like modulation 
since the correlation lengths become very small. Hav- 
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ing calculated the ground state in the corresponding S z - 
subsector we use equation ([h|) to deduce the improved 
set of {Si}, This step is repeated (typically 6 to 10 times) 
until the change of the ground state energy becomes suf- 
ficiently small, i.e. of the order of the truncation error. 

For the XY-Model we can compare the selfconsistent 
DMRG results for finite chains with these of the contin- 
ued fraction technique |l9[| in the thermodynamic limit. 
We find that the error due to finite size effects and due 
to the truncation of the Hilbert space is at most of the 
order of 10 -4 (see upper part of Fig.^). This accuracy is 
by far sufficient for the presented qualitatively analysis 
of the phase transition. 



IV. SUMMARY AND DISCUSSION 

In this work we considered modulated S = \ Heisen- 
berg chains with finite magnetization. Three classes 
of modulation were investigated: (i) fixed dimerization, 
(ii) fixed incommensurate sinusoidal modulation and (iii) 
adaptive incommensurate modulation. Our main interest 
was to investigate the dependence of the magnetization 
m on the applied magnetic field H. 

For scenario (i) we found a second order transition to 
finite magnetization by means of the finite size DMRG 
method in agreement with previous calculations. The in- 
crease of m just above the critical field H c is characterized 
by a square root behavior m ot yj H — H c . 

For scenario (ii) we showed that the system favors an 
incommensurability corresponding to the magnetization 
q = 2fcp = 27r(l/2 + m). A prominent first order tran- 
sition is found. This finding could be explained by com- 
puting the discontinuous dependence of the ground state 
energy E(m) on m at m = 0. The discontinuity is linked 
to the discontinuous jump of the root-mean-square of the 
local distortions on passing from dimerization (q — tt) to 
a long- wave length modulated dimerization (g w ?r). 

In scenario (iii) we determined iteratively the mod- 
ulation which minimizes the total energy including a 
quadratic elastic energy. Again we find a periodicity cor- 
responding to the magnetization q = 2n(\/2 + m). The 
modulation, however, corresponds for low magnetization 
rather to a soliton lattice. This means one has differently 
dimerized regions separated by domain walls. Each do- 
main walls carries one S — 1/2. We find a crossover from 
the solitonic picture at low magnetization to a sinusoidal 
modulation at higher magnetizations. The transition to 
finite magnetization is second order although the increase 
is exponentially steep. The inclusion of a positive disper- 
sion of the elastic energy alters the order of the transi- 
tion. It is first order then. An exponential attraction of 
the solitons was identified. 

As another interesting quantity we calculated the adia- 
batic gaps A + y_ corresponding to the increment (decre- 
ment) of the magnetization by unity. The independent 
determination of these gaps requires the complete knowl- 



edge of E(m). The calculation was also done for scenario 
(iii). It was shown that these gaps are finite in the adia- 
batic treatment. 

The second order phase transition in the com- 
mensurate case (|) is in agreement with the fact 
that measurements under applied field for instance on 
Cu2(C5Hi2N2)2Cl4 show no hysteresis effects ||. This 
substance is found to be an antiferromagnetic Heisen- 
berg ladder which is equivalent to a strongly dimerized 
quasi-one-dimensional Heisenberg chain. The magneti- 
zation increases continuously. The expected square root 
behavior near H c , however, was not observed. 

We do not know of a substance which can be described 
by pure sinusoidally modulated exchange couplings. The 
modulation in the incommensurate phase of CuGe03 is 
in fact very close to a sinusoidal modulation 10 llj]. Re- 
cently Lorenz et al. 32 measured the magnetic field 
dependence of the spontaneous strain e(H) in CuGeO-3 
which is in first approximation proportional to the elastic 
energy associated to the lattice distortion. It decreases 
very fast near H c and saturates approximately at 1/4 of 
the value in the dimerized phase for H « 22T. One can 
conclude that there is a crossover from a solitonic dis- 
tortion pattern for small magnetizations to a sinusoidal 
one for larger magnetizations. Our model allows - for 
parameter values reasonable for CuGe03 within a one- 
dimensional approach a = 0.35, K ss 18 (=> 5 = 0.014 
in the D phase J33-|3^] ) - to describe the above crossover 
quantitatively p2| . 

The feature so far not understood in CuGeOa is the 
first order phase transition D— >I. From our findings it is 
tempting to attribute this weak first order property to 
a positive dispersive elastic energy, i.e. a dip in uj(k) at 
k = 7T. Unfortunately, there is no experimental indica- 
tion for such a feature in the phonon spectra |36|. The 
spring constant K(k), however, in the adiabatic treat- 
ment is proportional to u>(k)/ g 2 (k) where uilk) is the 
phonon energy measured by Braden et al. ]3q] . The 
momentum-dependent spin-phonon coupling g^k) is not 
known presently and may account for the dispersive be- 
havior needed to explain the first order transition D— >I 
observed in CuGe03. At the present stage, we may also 
speculate that the neglected interchain couplings |37| are 
decisive for the order of the transition. From our present 
results we understand that the order of the transition is 
influenced by the microscopic details of the model. 

The finding of finite adiabatic gaps in the incommensu- 
rate phase should encourage experimental work to verify 
or to falsify this feature, for instance, in CuGeOs. 
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